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Abstract. In the framework of an SU(3) (axial)vector meson extended linear sigma model 
with additional constituent quarks and Polyakov loops, we investigate the effects of (axial)vector 
mesons on the chiral phase transition. The parameters of the Lagrangian are set at zero 
temperature and we use a hybrid approach where in the effective potential the constituent 
quarks are treated at one-loop level and all the mesons at tree-level. We have four order 
parameters, two scalar condensates and two Polyakov loop variables and their temperature and 
baryochemical potential dependence are determined from the corresponding field equations. We 
also investigate the changes of the tree-level scalar meson masses in the hot and dense medium. 


1. Introduction 

Forthcoming heavy ion experiments, such as the planned CBM experiment at FAIR will explore 
the QCD phase diagram in the high density and moderately high temperature region. This 
region is very interesting since it is believed that the critical endpoint (CEP) - which separates 
the crossover and first order phase transition regions along the phase boundary - if it exists, 
should be found there. 

Earlier experiments such as RHIC, SPS, or LHC investigated the low density, high 
temperature region. No sign of the chiral phase transition has been observed experimentally, 
which is not so surprising considering that the phase transition is of crossover type there, which 
is very hard to observe experimentally. Theoretically, if the transition is of first order - and the 
system is infinite - certain quantities have discontinuities crossing the phase boundary, which 
gives a much higher chance to observe it. Although, since in real physical experiments the 
system is always finite, there will not be any observable discontinuities, but just peaks. 

Due to the difficulties of the observation it is very important to gather as much theoretical 
information as possible from the phase boundary, in order to assist future experiments. 
According to that, we would like to investigate the chiral phase transition in the framework 
of a (axial)vector meson extended linear sigma model with additional constituent quarks and 
Polyakov loops. 

The paper is organized as follows. The model is briefly presented in Sec. [2] and its 
parametrization is discussed in Sec. [3j The four field equations which provide the temperature 
and chemical potential dependence of the order parameters are given in Sec. [4j Finally, the 
results are discussed in Sec. [5j where we also conclude. 


2. The Model 

In terms of matrix valued (pseudo)scalar and (axial)vector meson fields the Lagrangian is, 
C = Tr[(T> M $) t (D^)] - mg Tr^$) - Ai [Tr(& 4 >)] 2 - A 2 Tr($ f $) 2 


- \ Tr(T 2 „ + R 2 




/XI/. 


+ Tr 


^ + A)(L 2 +X ? 2 


+ Tr[XX($ + $+)] 


+ ci(det 4> + det $t) + i^(Tr{L^[L», L u ]} + Tr {R^[R^ XT]}) 

+ y Tr($t$) Tr(L 2 + X? 2 ) + h 2 Tr[(L M $) 2 + ($XX M ) 2 ] + 2X 3 TV(L M W$t) 

+ g^L„L v lJ x L v ) + Tr(R^R u R^R u )\ + 54 ['Tr (L^L V L U ) 

+ Tr (X^I^IT)] + Tr (L^XY) T4 (X^XT) + 5 6 [Tr(X^XX t ) Tr(L u L v ) 

+ Tr(X^) Tr^XT)] + #z 7/i XX^ - g F V ($ s + z 75 $ps) *, 


(1) 


where £>**$ = <9^4* - igi(LV<5> - <&R^) - ieA£[T 3 ,§], D w = & i L v - ieA£[T 3 ,L u ] - {d v - 
ieA u e [T 3 ,L^}}, R> w = d^R v - ieA%[T 3 ,R u ] - {d v R» - ieA u e [T 3 , R»}}, and XT 4 = <9^ - iG **. The 
held content of the Lagrangian is as follows, <f> is the scalar/pseudoscalar held, and R^ are 
the left and right handed vector helds, 'L = (u, d, s) T stands for the constituent quark helds, 
G^ is the gluon held, while XX is the external held. As usual, the nonstrange and strange scalar 
helds are shifted by their expectation values </>jv and (f>s (scalar condensates). A previous version 
of this model, without the constituent quarks and Polyakov loops and with a different anomaly 
term was soundly analyzed at zero temperature in [TJ. 


3. Determination of the parameters of the Lagrangian 

In Eq. 0 , not considering the unused 53,54,55, and ge , there are 14 unknown parameters, 
namely ?no - the bare (pseudo)scalar mass, Ai, and A 2 - the (pseudo)scalar self-couplings, ci - 
the 1 X^( 1 ) anomaly coupling, m 4 - the bare (axial)vector mass, h\, X 2 and h 3 - the (axial)vector- 
(pseudo)scalar couplings, 5s - the (axial)vector explicit break coupling, qbpr and 4>s - the scalar 
condensates, g F - the Yukawa coupling, and 51 and 5 2 - two (axial)vector couplings. Compared 
to the parametrization done in |T], in the present case we have two modifications: i) there are 
two additional equations for the constituent quark masses, m u /,i = g F (j)isr/2 and m s = g F 4> s /V 2 , 
for which we use the values m u /d = 330 MeV and m s = 500 MeV; ii) the inclusion of the 
fermion vacuum fluctuations modifies at T = 5 = 0 the masses and decay widths of the model 
from which the parameters were calculated in [1]. This is because in the present approach the 
masses used in the parametrization are the curvature masses which are the second derivatives of 
the grand canonical potential with respect to the meson helds at the minimum. The curvature 
masses consist of two parts: a usual tree-level part, which can be found in [TJ, and the vacuum 
and thermal contributions of the one-loop fermion potential (bosonic fluctuations are neglected). 
The later one can be found in mm with Polyakov loops, and in [ 3 ] without Polyakov loops. 
Scanning through the parameter space, curvature masses and decay widths are calculated and 
compared using a y 2 minimization method [5j to the experimental data taken from the PDG [ 6 ] 
and the above values of m u /d and m s , as described in [T], 

It is important to note, that in the scalar sector there are more physical particles than we can 
describe with one qq nonet, as in nature there are two ao, two Kq and five /o particles (see 0). 
Since in one scalar nonet there are one ao, one XLq and two fo’s, we have 2 • 2 ■ ( 2 ) =40 different 
possibilities for matching the scalar sector with physical particles. However, we consider here 
only cases where ao and XLq correspond to the ao(980) and XLq (800) physical particles. 

We apply two different parametrization scenarios here. In the first one we do not fit the very 
uncertain isoscalars (/q, /$)■, consequently mo and Ai always appear in the same combination 
C\ = mg + Ai {(p 2 N + (j )|) in all the expressions, thus we can not determine them separately. A 




Table 1. Parameters determined by x' 2 minimization in the two cases 


Parameter 

Value (16) 

Value (2) 

Parameter 

Value (16) 

Value(2) 

c^N [GeV] 

0.1359 

0.1333 

h 2 

4.8765 

2.7065 

<t>s [GeV] 

0.1400 

0.13823 

h 3 

4.6523 

3.7935 

m,Q [GeV 2 ] 

-0.0103 

0.0394 

5s [GeV 2 ] 

0.1114 

0.1178 

m 2 [GeV 2 ] 

0.5600 

0.5508 

ci [GeV] 

1.5293 

1.600 

Ai 

0 (undetermined) 

-1.2200 

9i 

5.5737 

5.5761 

A 2 

23.0638 

23.7957 

92 

2.1263 

1.3889 

hi 

0 (undetermined) 

2.6007 

9F 

4.3650 

4.4217 


similar combination C 2 = rn\ + ( <f>% + appears in the vector sector for mi and h\ [T]. 
Practically, we choose Ai = hi = 0 in this scenario. Even if we do not fit the isoscalars, we 
still have the possibility to get parametrizations with different values of mji, which has a huge 
effect, as will be seen immediately, on the finite T/ps behavior. We consider two subcases, 
labeled by ’la’ and T6’: one with a high and one with a low value of rrij-L, that is 1326 MeV 
and 402 MeV, respectively. In the second scenario, labeled as case ’2’, we fit the isoscalars to 
/o(500) and /o(1370). The parameter values for case T6’ and case ’2’ are given in Table [lj 


4. T / ps dependence of the order parameters and curvature masses 

In medium, the values of the order parameters - which are the two scalar condensates 4>n 
and 4>s and the two Polyakov loop variables <J> and $ - change with the temperature/chemical 
potential. The two scalar condensates encodes the effect of both spontaneous and explicit 
symmetry breakings, while the Polyakov loop variables mimic some properties of the quark 
confinement, which naturally emerge in mean field approximation, if one calculates free fermion 
grand canonical potential on a constant gluon background (for more details see [7j). 

In order to determine the T / ps dependence of the order parameters and curvature masses, we 
use four coupled stationarity equations (field equations), which require the vanishing of the first 
derivatives of the grand canonical potential with respect to the order parameters. We apply here 
a hybrid approach, where we only consider vacuum and thermal fluctuations for the fermions 
and not for the bosons. In this case the equations are given by 
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where t/( < f > , «4>) is the Polyakov loop potential, for which we used a polynomial form with 
coefficient taken from jS], and 


g+(p) = 1 + 3 (l* + $ e -^ +(p) ) e -^ + (P) + e —3 0E+{p)^ 
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Figure 1. Temperature dependence of the 
order parameters with the parametrization 
’la’ (to jl = 1326 MeV). 



T [MeV] 

Figure 2. Temperature dependence of the 
order parameters with the parametrization 
'lb' ( m f L = 402 MeV). 


E q(p) = E q{p) -F /W3, E u/d (p) = \jp 2 + m 2 u/d , E s (p ) = v / p 2 + m 2 , 
(qq) T = -4 m q J (l - f$(E q (p)) - f£(E q (p))) , 


with the modified Fermi - Dirac distribution functions 
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5. Results and Conclusion 

By solving the system of Eqs. (©-ft the temperature and baryochemical potential dependence 
of the order parameters and the curvature masses can be determined. In Fig. [I] the order 
parameters can be seen with the parametrization ’la’ at pb = 0. Here the pseudocritical 
temperature (T c ) is around 550 MeV, which is much higher than the continuum lattice result 
[9], which is T c = 151 MeV. In Fig. [2] the order parameters are shown with the parametrization 
’15’, where to^l = 402 MeV. In this case T c is between 150 — 200 MeV, which is in the range of 
the lattice results. 

It is a common belief that there is a critical endpoint (CEP) along the phase boundary in 
the T — pb plane. For this to happen, the phase transition should be of first order as a function 
of pb along the T = 0 axis. However, as can be seen in Fig. [3] (case ’16’) the phase transition 
at T = 0 is of crossover type, while in Fig. [4] with parametrization ’2’ the transition is of first 
order. In this case the value of the pseudocritical temperature T c at ps = 0 is close to that of 
case 'lb'. The temperature dependence of the pseudo(scalar) curvature masses can be seen in 
Fig. [5] and Fig. [6] for the parametrization of case ’2’. 

In conclusion we can say that the parametrization which fulfills the two physical requirements 
of having a pseudocritical temperature T c ~ 150 MeV at pb = 0 and a first order transition in pb 
at T = 0 is realized in case ’2’, that is when /o(500) and /o(1370) are used for parametrization. 
Similar results were found in Ell without (axial)vector mesons. 














Figure 3. fj ,b dependence of the condensates 
with the parametrization ’16’. 



Figure 5. Temperature dependence of the 7r, 
r/, ao and a) particle masses. 



Figure 4 . hb dependence of the condensates 
with the parametrization ’2’. 



Figure 6. Temperature dependence of the K , 
r/, K* and particle masses. 
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